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i Unit 30480685, "Rotor-Bearing Dynamics Design." 


The work reported herein was performed during the period 15 September _ I: 
1977 to 15 June 1978, under the direction of John B. Schrand (AFWAL/ 
POSL and Dr. James F, Dill (AFWAL/POSL), Project Engineers. The report 
was released by Shaker Research Corporation in January 1979, 


TABLE OF CONTENTS 


Section INTRODUCTION 
Section 
Section ITI 


3.1 
3.2 
3.3 


Section 


4.1 
4.2 
4.3 


Section 


5.1 
5.2 
5.3 Reynolds Condition 


Conclusions 


REFERENCES 


LIST OF ILLUSTRATIONS 


Bearing Coordinate System 


Single Pad Schematic 


Principal Stiffnesses vs. Sommerfeld Number for Plain 
Journal Bearing 


Cross-Coupled Stiffnesses vs. Sommerfeld Number for 
Plain Journal Bearing 


Principal Damping vs. Sommerfeld Number for Plain 
Journal Bearing 


Cross-Coupled Damping vs. Sommerfeld Number for Plain 
Journal Bearing 


Stability Threshold vs. Sommerfeld Number for Plain 
Journal Bearing 


Principal Stiffnesses vs. Sommerfeld Number for Plain 
Journal Bearing 


Cross-Coupled Stiffnesses vs. Sommerfeld Number for Plain 
Journal Bearing 


Principal Damping vs. Sommerfeld Number for Plain 
Journal Bearing 


Cross-Coupled Damping vs. Sommerfeld Number for Plain 
Journal Bearing 


Table 


10 


11 


12 


13 


LIST OF TABLES 


Plain Journal Bearing Dimensionless Characteristics __ 


Two Lobe Bearing Dimensionless Characteristics 

L/D = 0.5, m= 0.0, x = 160, Load on Pad 
Two Lobe Bearing Dimensionless Characteristics 

L/D = 0.5, m= 0.25, x = 160, Load on Pad 
Two Lobe Bearing Dimensionless Characteristics 

L/D = 0.5, m = 0.5, x = 160°, Load on Pad 


Three Lobe Bearing Dimensionless Characteristics 

L/D = 0.5, m= 0.0, x = 100°, Load Between Pads __ 
Three Lobe Bearing Dimensionless Characteristics 

L/D = 0.5, m = 0.250, x = 100, Load Between Pads __ 
Three Lobe Bearing Dimensionless Characteristics 

L/D = 0.5, m = 0.50, x = 100°, Load Between Pads 
Dynamic Data of Fixed Pad (L/D = 0.25, Offset = .500, 

x = 80°, Load on Pad Center) 
Static Data of Fixed Pad (L/D = 0.25, Offset = .500, 

x = 80°, Load on Pad Center) 
Four Shoe Tilting Pad Bearing (L/D = 0.25, m = 0, Off- 
set = .500, x = 80 , Load Between Pads) Dimension- 

less Synchronous Dynamic Data 

Four Shoe Tilting Pad Bearing (L/D = 0.25, Offset = 
500, m = 0.2, x = 80°, Load Between Pads) Dimen- 
sionless Synchronous Dynamic Data 

Four Shoe Tilting Pad Bearing (L/D = 0.25, Offset = 
-500, m = 0.3, x = 80 , Load Between Pads) Dimen- 
sionless Synchronous Dynamic Data 

Four Shoe Tilting Pad Bearing (L/D = 0.25, Offset = 
-500, m = 0.5, x = 80°, Load Between Pads) Dimen- 
stonless Synchronous Dynamic Data 


ee ee ee ew ee ee ee ee 


20 


21 


22 


23 


24 


25 


26 


28 


29 


30 


31 


32 


33 


14 


15 


16 


17 


18 


19 


20 


Ses oe eee 
ve, 


21 


22 


23 


24 


25 


LIST OF TABLES 


(Continued) 


Dynamic Data of Fixed Pad (L/D = 0.50, Offset = 
-500, x = 80° Load on Pad Center) 


Static Data en Pad (L/D = 0.5, Offset = 
-500, y = 80°, Load on Pad Center) 


Four Shoe Tilting Pad Bearing (L/D = 0.5, Offset = 
-500, m = 0.0, x = 80° , Load Between Pads) Dimen- 
sionless Synchronous Dynamic Data 


Four Shoe Tilting Pad Bearing (L/D = 0.5, Offset = 
-500, m = 0.2, x = 80, Load Between Pads) Dimen- 
sionless Synchronous Dynamic Data 


Four Shoe Tilting Pad Bearing (L/D = 0.5, Offset = 
.500, m = 0.3, x = 80, Load Between Pads) Dimen- 
sionless Synchronous Dynamic Data 


Four Shoe Tilting Pad Bearing (L/D = 0.5, Offset = 
-500, m = 0.5, x = 80 , Load Between Pads) Dimen- 
sionless Synchronous Dynamic Data 


Dynamic Data of Fixed Pad (L/D = 1.00, Offset = .500 
x = 80, Load on Pad Center) 


se Data of Fixed Pad (L/D = 1.0, Offset = .500, 
= 80°, Load on Pad Center) 


Four Shoe Tilting Pad Bearing (L/D = 1.0, Offset = 
500, m = 0.0, y = 80 , Load Between Pads) Dimen- 
sionless Synchronous Dynamic Data 


Four Shoe Tilting Pad Bearing (L/D = 1.0, Offset = 
-900, m= 0.2, x = 80° ,» Load Between Pads) Dimen- 
sionless Synchronous Dynamic Data 


Four Shoe Tilting Pad Bearing (L/D = 1.0, Offset = 
-500, m = 0.3, xy = 80 , Load Between Pads) Dimen- 
sionless Synchronous Dynamic Data 


Four Shoe Tilting Pad Bearing (L/D = 1.0, Offset = 
-500, m= 0.5, y = 80 , Load Between Pads) Dimen- 
sionless Synchronous Dynamic Data 


ee ee re ee ee ee 


ee ee ae we ee ee ee 


er ww nw ow ee we ee ow ee = 


35 


36 


37 


38 


39 


40 


41 


42 


43 


44 


45 


LIST OF TABLES 


(Continued) 
Page 
26 gree: Data of Fixed Pad (L/D = 0.25, Offset = .500, 
= 55°, Load on Pad Centers ooh ist i he 46 
27 ee Data of Fixed Pad (L/D = 0.25, Offset = .500, 
= 55°, Load on Pad COREE ce ao ete ie de al cg 47 
28 Five Shoe Tilting Pad Bearing (L/D = 0.25, Offset = 
-500, m = 0.0, x = 55°, Load on Pad) Dimensionless 
Synchronous Dynamic Data 48 
29 Five Shoe Tilting Pad Bearing (L/D = 0.25, Offset = 
500, m = 0.200, x = 55 , Load on Pad) Dimension-~ 
less Synchronous Dynamic Data_ 49 
30 Five Shoe Tilting Pad Bearing (L/D = .25, Offset = 
-500, m = 0.300, x = 55°, Load on Pad) Dimension~ 
less Synchronous Dynamic Data 50 
31 Five Shoe Tilting Pad Bearing (L/D = 0.25, Offset = 
-500, m = 0.50, xy = 55°, Load on Pad) Dimension- 
less Synchronous Dynamic Data 51 
32 Dynamic Data of Fixed Pad (L/D = .50, Offset = .500, 
= 55°, Load on Pad Center) 52 
33 Static Data of Fixed Pad (L/D = 0.5, Offset = .500, 
= 55°, Load on Pad Center) 53 
34 Five Shoe Tilting Pad Bearing (L/D = 0.5, m= 0.0, 
= 55°, Load on Pad) Dimensionless Synchronous 
Dynamic Data’... \2 si oe ee ee ee 54 
35 Five Shoe Tilting Pad Bearing (L/D = 0.5, Offset = 
-500, m= 0.2, x = 55°, Load on Pad) Dimensionless 
Synchronous Dynamic Data- 55 
36 Five Shoe Tilting Pad Bearing (L/D = 0.5, Offset = 
-500, m = 0.3, y=55 , Load on Pad) Dimensionless 
Synchronous Dynamic Data 56 
37 Five Shoe Tilting Pad Bearing (L/D = 0.50, Offset = 
-500, m= 0.5, y = 55°, Load on Pad) Dimensionless 
Synchronous Dynamic Data 57 


ix 


ae ee ee oe we ee ow em ne ee eee 


LIST OF TABLES 


(Continued) 


38 Dynamic Data of Fixed Pad (L/D = 1.00, Offset = 
-500, x = 55°, Load on Pad Center) 
39 Static Data of Fixed Pad (L/D = 1.0, Offset = 
-500, x = 55', Load on Pad Center) 
40 Five Shoe Tilting Pad Bearing (L/D = 1.0, Offset = 
-500, m = 0.0, x = 55 , Load on Pad) Dimensionless 
Synchronous Dynamic Data 
41 Five Shoe Tilting Pad Bearing (L/D = 1.0, Offset = 
-500, m = 0.2, x = 55 , Load on Pad) Dimensionless 
Synchronous Dynamic Data 


we ew ee ee ewe ee we we ee oe we 


42 Five Shoe Tilting Pad Bearing (L/D = 1.0, Offset = 
-500, m = 0.3, x = 55 , Load on Pad) Dimensionless 
Synchronous Dynamic Data 
' 
: 43 Five Shoe Tilting Pad Bearing (L/D = 1.0, Offset = 
iG -500, m = 0.5, x = 55, Load on Pad) Dimensionless 
‘ Synchronous Dynamic Data 
: 
INDEX FOR TILTING PAD BEARING DATA TABLES 
’ 
AP 


Symbol 


LIST OF SYMBOLS 


A. Dimensional Quantities 


Definition Units 
Damping coefficient lb-sec/in 
Clearance (nominal or bearing radial in 


clearance for plain journal; pad or 
ground in radial clearance for multi- 
lobe and tilting pad bearings) 


Bearing clearance (radial clearance of in 
largest cylinder that will fit in bear- 
ing, preload m = (€ - G,)/c) 


Diameter of shaft in 
Eccentricity in 
Bearing eccentricity in 
Bearing film thickness in 
Pivot film thickness in 
Horsepower loss hp 
Stiffness Teta 
Length of bearing 3 
Mass of the rigid rotor lb-sec’ /in 
Rotational speed of shaft rev/min 
Rotational speed of shaft rev/sec 
Pressure lbf/in 
Inlet oil flow to bearing (sum of inlet in?/sec 


oil flows to all pads in multilobe and 
tilting pad bearings) 


Lost oil flow from bearing (sum of axial in?/sec 
oil flows of all pads in multilobe 
and tilting pad bearings) 


Radius of preload circle, tp =Cc-C¢ in 


Journal radius 


Pad radius of curvature in 


xi 


LIST OF SYMBOLS 


, 
(Continued) : 
Ry Radius from bearing center to pivot in 
Time sec f 
T Friction torque 1bf-in 
Load on bearing (vertically downward lbf. 
steady state load) 
x, Coordinate system for bearing in 
: z Axial coordinate in 
ra not Fixed pad coordinate system degrees : 
a 5 Pitch angle of pad degrees 
i 8 Angular coordinate (from Y-axis) degrees 
. 8 Angle from leading edge of pad to pad degrees 
= P pivot 
7 Pad arc length degrees 
.* 
- Angle from positive y-axis, degrees 
PS 1 ith pad pivot point measured in direction 
,e of rotation measured from positive y-axis 
f w Rotational speed rad/sec } 
Average viscosity Upeceee/ ine I 
q v Attitude angle measured with rotation degrees i 
5. from the static load vector (positive ; 
a X-axis) to the eccentricity vector i 
B. Dimensionless Quantities | 
| 
’ Symbol 
B Bew/WS, dimensionless damping coefficient | 
h,h h_/c_, dimensionless pivot film thick- { 
p’ pi p Pp th 
, ness, pivot film thickness for i- pad 
K Ke/WS, dimensionless stiffness coefficient 
= 22.3 2 pas ; 
M 4n MC No/uR LD, critical mass parameter | 
m 1 - (G/0) bearing preload factor 
n Number of pads 


LIST OF SYMBOLS 


(Continued) 


p -2 
K +5 
nn Ban 


Q/N,LDc, dimensionless oil flow parameter 
uNLD(R/c)“/ (60W), Sommerfeld number 
H (4) (12) (550)¢/2*N, “ LD”) 


Dimensionless horizontal, vertical hydrodynamic load 


Dimensionless horizontal, vertical hydrodynamic load 
for the i- pad 


6! % bearing offset factor 


e,/6, eccentricity ratio 
e,/C,> eccentricity ratio 


iB ne 


SECTION I 


INTRODUCTION 


This work is intended to supplement and in some cases supersede the 

report entitled "Rotor-Bearing Dynamics Design Technology, Part III: 

Design Handbook for Fluid Film Type Bearings," Technical Report 
AFAPL-TR-65-45, Part III, May 1965. This report discussed static and dynamic 
properties of plain cylindrical bearings, four axial groove bearings, 
elliptical bearings, partial arc bearings, and tilting pad bearings. This 
report was followed by another entitled "Rotor-Bearing Dynamics Design 
Technology, Part VII: The Three-Lobe Bearing and Floating Ring Bearings," 
Technical Report AFAPL-TR-65-45, Part VII, February 1968. As indicated by 
the title, this report extended the 1965 report by including three-lobe 


journal bearings and floating ring bearings. 


The subject of the present study is to augment the bearing data given 

in the two previous reports by including a significant amount of additional 
multilobe bearing data, including several preload values. Preload values 
of 0.0 (a three-axial groove bearing) and 0.25 are considered in addition 
to the previously considered values of 0.5. The present report also 


supersedes the bearing data in the previous report AFAPL-TR-65~45, Part III 


for the tilting pad data. Apparently some inaccuracies were present in the 


pad assembly computer program utilized to calculate the data for the original 


curves. This report presents a wide variety of tilting pad static and dynamic 


properties for calculational use. 


Comparisons of the data in this volume and a number of published works gives 
good correlation for cases where identical bearings have been considered by 
different authors. Generally, the errors are on the order of magnitude of 


10% or less. There will always be minor differences from one set of 


numerical calculation to another due to the numerical “technique used (finite 


d 
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elements versus finite differences, the iteration procedure used to find 
equilibrium position, etc.) such that errors of the 5 to 10% order of magnitude 
are reasonable. It has also been found by many authors that the correlation 

of theoretical with experimental results is usually much worse than the 5 to 


10% range. Thus, errors of the 5 to 10% range are not critical. 


The data given in this report includes the equilibrium position, bearing 
coefficients, horsepower loss, and oil flow of the following bearing types: 


» plain journal 


multiiobe | 
tilting pad 

A finite element solution of the Reynolds equation is used to obtain the pres- 

sure over the finite bearing surface. Pressures that are lower than the 

cavitation pressure are set equal to the cavitation pressure. The effects 

of fluid streamers forming in the cavitation region are included. The 

lubricant is assumed to be isoviscous and laminar. The pressure is then 

integrated over the bearing to obtain the load capacity. For a given load 

the equilibrium eccentricity, and attitude angle are determined by the New- 

ton Raphson iterative procedure. Linear stiffness and damping coefficients 

are obtained from small perturbations of displacement and velocity about the 


equilibrium position. The effects of misalignment of the shaft in the bear- 


ing are not considered. 
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SECTION II 


LUBRICATION THEORY 


The coordinate system (fixed in space) adopted in this report is shown in 


Fig. 1. Face the bearing each time so that the shaft rotates counterclockwise. 


The angular coordinate, 9, is measured positive in the counterclockwise direc- 


tion starting always from the positive y-axis. x and y are the vertical 


and horizontal coordinates of the geometrical center of the shaft (centerline). 


The form of the Reynolds equation solved is 


sy 


3 3 \ 


1 a {hap |, ap {n° ae | Gy ah, 12n dh 
RZ 20 90 | z az 30 at (1) 


All of the standard Reynolds equation assumptions hold. 


1. The average clearance is small compared to the length in either 


the x or z directions. 


2. Only laminar flow occurs in the bearing. The Reynolds number 
is less than 2,000. 


3. No-slip boundary conditions apply at fluid-solid interfaces. 


4. All body or other external forces are neglected in the fluid 
film. 


5. The fluid is Newtonian, incompressible, and isoviscous. 


The boundary condition used is the half Sommerfeld condition. Pressures be- 
low the cavitation pressure, (zero psi) are set equal to the cavitation pres- 


sure and it is assumed that fluid streamers form in the cavitated region 


(See Fig. 2). The determination of which areas of the bearing are cavitated 


is determined element by element. 
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Figure 1 Bearing Coordinate System 


+ 
Streamers in Cavitated Region (Top View). 
Figure 2 Boundary Conditions for Plain Journal Bearing 


Symmetry im the z-direction is assumed and the bearing split in half along 
its length (along z = 0) as shown in Fig. 3. Then, the hydrodynamic pres- 


sures are solved for only half of the bearing, saving computer time. 


For a plain journal bearing, the bearing is split at the point of maximum 
film thickness. The hydrodynamic pressure is assumed zero at hoax’ Accord- 
ing to Pinkus and Sternlicht, a full journal bearing is not very sensitive 

as to where the hydrodynamic pressure is assumed zero as long as it is 

near hae They state that even shifts of the order of 20° introduce little 
error. They observe that the performance of a 180° bearing is seen to differ 
little from a full 360° bearing. A full discussion of the boundary conditions 
is given in Section 5. 


For pad bearings, the hydrodynamic pressure is again assumed zero at 
z= + L/2. Symmetry in the z-direction is also assumed and the hydro- 


dynamic pressure is assumed zero at both ends of the arc. 


naax or leading edge haax or trailing edge 


Figure 3 Boundary Conditions 


SECTION IIT 


BEARING TYPES 


Fluid films are found in many rotating machines either as bearings, seals, 
or squeeze film dampers. Certain standard geometries such as those dis- 
cussed here have evolved to satisfy various requirements. Because bearings 
or seals are much more easily changed than other rotor parts, clearances, 
preload factors, or bearing types are often varied as vibration problems 
occur. It is desirable to be able to analyze the effects of those param- 


eters of all of these bearing types. 


3.1 Plain Journal Bearing 


As shown in Fig. 4, the plain journal bearing is composed of a cylindrical 
shaft of radius R, called a journal, rotating with angular velocity, w, 
counterclockwise about its axis in a cylindrical bushing of radius R + C 

and length L. The center of the journal is labeled Oy. and the center of the 
bushing is labeled Oo, as shown in Fig. 4. If the bearing is operating 
under steady state conditions, the journal center remains at a constant ec- 
centricity e and attitude angle y for a given load W acting on the shaft. 
The case of the plain journal bearing is an important one because this 
geometry is the same as that in a plain seal and also in a squeeze film 
damper (where w= 0). For this reason the analysis remains important al- 
though the use of plain journal bearings is declining due to stability 


problems. 


3.2 Multilobe Bearing 


The multilobe bearing is different from the plain journal bearing in that 
the pad of radius, R+C, has the center of curvature, On located a 


small distance T from the bushing center. The preload factor m equals 


Line of Centers 
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Figure 4 Plain Journal Bearing 
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to the pad surface at point P gives the 


r/¢. Extending the line 09 


location of minimum film oetecaae Che This corresponds to the maximum 

shaft radius R + C. which will fit in the bearing. The point P also 
determines the end of the converging portion of the pad and the offset 
factor « = a/x. Both the preload factor and the offset factor have 

values between zero and one. The offset factor is often chosen as one-half 


to permit running in both directions, Fig. 5 shows the geometry. 


Multilobe bearings offer the user an increase in stability in some situa- 
tions by the increased pad convergence loading up the bearing. Because of 
the many possible combinations of number of pads, rotation of bearing, 

clearance, preload factor, and offset factor, a large range of load carry- 


ing capacity and bearing coefficients is possible. 


A special case occurs when the preload is zero. This is often called an 
axial groove bearing as it may be constructed by cutting axial grooves in 
a plain journal bearing. It is often used as an inexpensive replacement 
for plain journal bearings in an attempt to reduce instability problems in 


a machine. 


3.3 Tilting Pad Bearing 


The tilting pad bearing shown in Fig. 6, differs from the multilobe bear- 
ing in that each pad is pivoted behind point P. Three or more pads are 
usually used. Each pad may have as many as four degrees of freedom corres- 
ponding to movement of the pivot point in or out due to a spring loading and 
three rotations. Of the three rotations, only pitch in the x,y plane is con-~ 
sidered here while roll and yaw are neglected. This corresponds to as- 
suming that the shoe follows the shaft as it precesses. In oil film bear- 
ings this approximation is usually a good one due to the high damping proper- 


ties. The approximation is not valid for gas bearings as pad inertia is 


important. 


Preload 
Circle, 
Radius ae 


Figure 5 Multilobe Bearing Geometry 


Largest Shaft 


Clearance, 
Radius R + C. 


R,,= RADIUS TO PIVOT 
R,= PAD RADIUS OF CURVATURE 
R= Rp WHEN my, = 0.0 
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FIGURE 6 TILT PAD BEARING SCHEMATIC 
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Tilting pad bearings are often used because they are more stable than other 
types of bearings. Because the pads are able to follow the shaft motion, 
there is little cross-coupzed stiffness and damping which is what drives 
fixed pad bearings unstable above certain speeds. The hydrodynamic load 
vector developed by each pad under steady state operating conditions must 


pass through the line between the pivot point and the journal center O.- 


For liquid lubricated bearings, the pad inertia may be neglected [Lund, 
1964]. With symmetry about the x-axis, the stiffness and damping coefficients 


for the full tilt pad bearing in dimensionless form are [Lund, 1964] 


n 

Ky (R’,,cos"4), (2) 
n 

ae > (RK, ,sin’ 4), (a) 
n 

Boy “ (B,,cos"4) , (4) 
n 

Be: “b, (Bein), (5) 

K =K =B  =B =0 (6) 


where n = number of pads 


, = angle from positive y-axis measured counterclockwise to the 


pivot point of the qe? pad (Fig. 6) 


and 


-_- 


BY. =B (PK ,, + @B ,,) B (qk 
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(8) 


with 


> 
nn q 
K + B 
nn nn ‘ 
i 
a i 
an : 
q = = (10) : 
x2 + 32 ‘ 
nn nn et 
s 
i 
f 
“+ The first subscript in the K and B terms refers to the force direction while : 
an the second refers to the direction of displacement or velocity perturbation. i 


The quantities on the right-hand side of Eqs. (7) through (10) 

7; [eens cette: cee 
: (Kee Kent She 

coefficients. 


Koy Bons, B the di i fixed pad 
ain? Bees Bey? Bue? =) are e dimensionless fix P 


These coefficients are calculated by loading up a single 


hil RAINE ATONE gO, HY WREIO EERE 


fixed pad, obtaining an equilibrium position and perturbating about the 
e equilibrium position to determine the fixed pad coefficients (see Fig. 7). 
e Then the pads dynamic coefficients, Khe and Bhe> are calculated using Eqs. (7) 
a through (10). The other 6 dynamic coefficients are zero because pad inertia is f 
 . neglected. These 2 dynamic coefficients along with the fixed pad Sommerfeld } 
1 number S are stored as a function of the pads dimensionless pivot film thick- \ 
' ness, hy where 


Sohne SAL ME 


h =h/C 
Pp Pp 
and 


5 = pad radial clearance 


i = pivot film thickness 


* ae 
Each B coefficient appearing within parentheses on the right-hand sides of 


Eqs. (7) and (8) and in Eqs. (9) and (10) should be multiplied by the fre- 
quency ratio v/w for nonsynchronous vibrations. 


Sewers ONE: Gate wee deel eS ak pee et WS Ate eo Wo eee eee 


Single Pad Schematic 
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Figure 7 


eae 
i 


To determine the bearing characteristics of the full tilt pad bearing, 


a bearing eccentricity ratio, ¢, is chosen where 


b 
*b 


e = c. (11) 
and 
Cy = tilt pad bearing assembled radial clearance in line with 
a pivot 
e, = tilt pad bearing eccentricity (Fig. 6) 


The bearing clearance C_, can be defined by the largest shaft ‘shaft 


b 
radius = R + Cy) which can be placed in the bearing after it is assembled. 


For a bearing with zero preload, the shaft with radius R + C, would be in 


b 
contact with the pad at all points. For a bearing with nonzero preload, 


shaft-pad contact will occur only along the pivot line. 


Secondly, the dimensionless pivot film thickness is calculated for the 
: P. Next Kea? Beri and S; are determined for the corresponding 
h_ value uSing a 4 point Lagrangian interpolation. This is repeated for 
eath pad. Eqs. (2) through (5) are then used to determine the dimen- 
sionless stiffness and damping coefficients for the full tilt pad bear- 


ing. Now, the dimensionless load must be calculated for each pad 


= -1 
= 2 
Wy S, cos o, (12) 
W@, = -s;tsin 6 (13) 
xi i i 
where 
5, = Sommerfeld number for hag pad 
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For the full bearing 


n 
fo) a (4) 
Y 21 Yi 
n 
Rnd a (5) 
i=l i 


The Sommerfeld number for the assembled tilt pad bearing is 


_as-l 
S=W, (16) 
and 
W = 0 
y 
The top pads of a zero preloaded tilt pad bearing have dimensionless pivot 4 


film thicknesses, hy » greater than 1.0 for the case of a horizontal shaft 
with a downward load. However, these pads still contribute to the bearing 


damping. That is, even if hy >1.0, B is nonzero and positive and contributes 


= g&. a 
to By and Beg in Eqs. (4) and (5). On*the other hand, Kee approaches 
zero as hy approaches 1.0. Thus, it is set to zero for ho 1.0. 

i i 


The assumption of symmetry about the x-axis holds for all load on pad (load 
on pivot) and load between pad (load between pivot) bearings. Also, the 
transformation from fixed pad data to dynamic pad data, Eqs. (7) through 
(0), may be viewed as an adjustment in the fixed pad data to allow the 

pad to pitch so that the load, W will pass through the pivot point. 


BEARING DATA 


| 
| 
SECTION IV 
| 
| 


Bearing data has been calculated for plain journal, multilobe, and tilting 
pad bearings. The data includes the following: 


a) Eccentricity, € 


tt 


€ = (plain journal) 


2 


4 ey = = (multilobe and tilt pad) 
b 
¢ b) Attitude angle, y 


c) Sommerfeld number, S 


F d) Torque, T 


4(12) (550)c 


P on N,uLD? 


T= 


e) Oil flow, Qreq? Qiost 


Q= pes. ee 
N LDC 


(Q 


Inlet flow to bearing 
req 


Qiost™ Axtal flow lost to side of bearing) 


f) Stiffness, Ky? Key 


K=K 


we 
WS 
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g) Damping, Bex? B ,B_,B 


eo. g Co 
B= B WS 
4.1 Plain Journal Bearings 


A plain journal bearing with L/D = 0.5 is considered in this section. 


The dimensionless characteristics are given in Table l. 


The value L/D = 0.5 is very common for journal bearings. Usually the 


range is from L/D = 0.4 to L/D 0.75 with not much variation within this 


range. Thus L/D = 0.5 represents most of the bearings in common use. 


4.2 Multilobe Bearings 


Multilobe bearings with L/D = 0.5 are presented in Tables 2 through 7. 
Both two-lobe and three-lobe bearings with preloads of 0.0, 0.25, 0.50 are 


evaluated. For the two-lobe bearing, the pad arc length is 160° and the 
load vector is through the center of the bottom pad. In the three-lobe 
case, the pad arc length is 100° and the bearing is rotated so that the load 
vector is directed between pads. Thus large loads are taken by two pads 
rather than just one. In all cases with non-zero preload, the offset factor 
(location of minimum film thickness, if the shaft is centered, expressed as 


a fraction of pad arc length) is 0.50. 


4.3 Tilting Pad Bearings 


Two tilting pad bearing cases are evaluated: four-pad with the load 
between the pads and five-pad with the load on the bottom pad. In each 
case, the L/D ratios are 0.25, 0.50, 1.00 and the preloads are 0.0, 

0.2, 0.3, 0.5. In all cases the pads are centrally pivoted. Page 27 con- 


tains an index for various data tables of tilting pad bearings. 
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SECTION V 


BOUNDARY CONDITIONS 


This section discusses the effects of the cavitation condition on bearing 
static and dynamic properties. An important point concerns the location 

of the boundary conditions during the calculation of dynamic bearing prop- 
erties. Specifically -- are they fixed at the same position obtained in 
Steady state operation or are the positions recalculated for each new co- 
efficient. The plain journal bearing is used as an illustration of the 
various results. It will be seen from the results that half Sommerfeld 
conditions with recalculated (floating) boundary locations during dynamic 
analysis gives good agreement with the Reynolds condition using fixed bound- 


ary locations in dynamic analysis. 


5.1 Types of Conditions 


A number of different boundary or cavitation conditions are currently 
used in bearing analysis. These may be generally grouped under the 


following headings: 


Full Sommerfeld 
Half Sommerfeld Boundary or Cavitation Conditions 


Reynolds 


The full Sommerfeld condition involves neglecting the effects of cavitation 
within the bearing. The half Sommerfeld condition consists of solving 
Reynolds equation (described in Section 2) and then setting the resulting 
negative pressures to ambient (zero) pressure. Finally Reynolds condition 
requires that both the pressure and the slope of the pressure normal to the 


boundary is zero at the trailing edge boundary. 


The effects of cavitation are very important in bearings so that full 
Sommerfeld condition is not considered in this report. Only under very 


high ambient pressure such as in a nuclear water pump does this condition 
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apply. In a plain journal bearing the full Sommerfeld condition gives an 


unstable bearing. 


In evaluating the half Sommerfeld condition, the steady state running condi- : 
tions are determined by assuming a vertically downward load and then iterating 
until the equilibrium angle is obtained. The location of the leading edge 
boundary is at the maximum film thickness in a plain journal bearing or 
diverging-converging pad while it is at the leading edge of a converging- 
diverging pad. The trailing edge boundary is at the minimum film thicness in 

a plain journal bearing or converging-diverging pad but at the trailing edge of 

a diverging-converging pad. During the dynamic analysis, the diaplacement 

and velocity perturbati zs to determine the dynamic coefficients produce small 
changes in film thickness, The following analysis of a plain journal bearing 
shows that these small changes in film thickness do not produce small (neglectable) 
effects in the half Sommerfeld condition dynamic bearing coefficients. If the 
conditions are kept at the steady state position — both leading and trailing 
edges — the conditions are called fixed-fixed in this report. If they are 
recalculated for each dynamic coefficient — both leading and trailing edges — 


the conditions are called floating-floating in this report. 


Reynolds condition is used to obtain the equilibrium position in a manner 
similar to that in the half Sommerfeld condition except that the trailing 
edge location is iterated over the pad surface until both the pressure and its 
derivative are zero. In the case of a purely converging pad, this procedure 
is ignored. Usually the effect of the boundary curvature is small so the 
angle at which Reynolds condition is applied may be considered constant. 

The boundary conditions may be either considered as fixed-fixed during 


dynamic calculations or as floating-floating. 


At the leading edge oy and the trailing edge 6 of the positive pressure 


region, the boundary conditions for the half Sommerfeld condition are 
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P=0 , 9 
P=0O , 8 


Le 


8s 


while for the Reynolds condition 


P=0 cy alae ae 
ap - 
Boe Oe oe 


These conditions refer to plain journal bearings but similar results can be 
applied to other bearings. The last equation indicates that the derivative 


of P normal to the trailing edge curve is zero. It may also be written as 


oP op 
76 Bo aces te 


where ng and n, are the direction cosines of the normal trailing edge 


curve. This is zero only if 


aP _ 
SG a? Oe 
oP = 

ee ea ee 


An additional subscript has been added to the leading and trailing edge 
angles to denote half Sommerfeld (s) and Reynolds (r) conditions. The 
half Sommerfeld condition always leads to a positive pressure region of 
n radians and a negative pressure region of 7 radians so that Use may be 


taken either as 


ee &,t 27, and P > 0 


The Reynolds boundary condition calculates a value for 63, such that 
T < Gi - Soe < 2m and neglects the pressure in the remainder of the 


bearing. 


An alternative method of calculating the dynamic coefficients to the 
numerical perturbations employed in this report is the direct perturba- 
tion of the differential equation [tund, 1966], The film thickness is 
analytically perturbed and a separate differential equation is solved 
(numerically) for the dynamic coefficients. In applying the Reynolds con- 
dition, it can be shown that the derivative of the perturbed pressure is 
zero at the same location as the steady state pressure (to first order). 


Thus the fixed-fixed conditions are normally used with Reynolds condition. 


The direct perturbation method assumes that the amplitudes of the first 


order perturbation pressures can be written as 


= + P + + P.Ax + P.AY 
P Po OX PyAy AX PLAY 


Pressure Under Static Equilibrium Conditions 
Pressure Perturbation to Ax Displacement 
Pressure Perturbation to Ay Displacement 
Pressure Perturbation to 4x Velocity 


Pressure Perturbation to Ay Velocity 


Applying the perturbation technique to Reynolds equation yields five equations 
which must each be solved for the pressure [Lund, 1966], (Lund, 1978]. The 
boundary conditions at the trailing edge (the other boundaries are discussed 
in the references and are the same for half Sommerfeld and Reynolds condi- 


tions) are subject to a first order expansion also 


- a) sor(82) ae « 
P(6,z) = P(8 320) + 35) aoe) Az 0 (18) 


Rewriting Eq. (17) as 


P® Py + AP 


Eq. (18) 
P(6,z) = Po (8qs 2) + AP(6,5 Z) 
oP oP 
0 0 
+ (zl, Ae + Fee 4z = 0 (19) 


It is apparent from this equation that for Reynolds condition the pres- 


sure gradients vanish 


(32°) = a =zQ@ 
a6 lo dz 0 


and the perturbation pressure must vanish at the same trailing edge as 

used in the static equilibrium state. However, for the half Sommerfeld 
condition, the derivative of P with respect to 6 is always large even 

if OP ae is small so Eq. (19) becomes 


aPy | oP 
AP(8); Zo) + 20 by AG + Pe Az = 0 (20) 
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(since P. = 0). The three terms of Eq. (20) are first order effects of the 


0 
perturbation and therefore must all be retained. Therefore, for the half 


Sommerfeld solution, the boundary conditions for the perturbed pressure may 
not be zero at the trailing edge. An equivalent numerical result is achieved 
if the trailing edge boundary be allowed to float during calculation 


of the dynamic coefficients. 


5.2 Half Sommerfeld Condition 


The half Sommerfeld condition was used to determine stiffness and damping 
coefficients for a plain journal bearing. Both fixed-fixed and floating- 
floating conditions were evaluated and compared to Reynolds fixed-fixed 
condition [Lund, 1966] . All of the stiffness and damping coefficients in 
this section were calculated using a series solution to Reynolds equation 
[Li, Allaire, Barrett, 1977] with a large number of series terms included 
so that no numerical errors are involved there. The effect of numerical 
displacement and velocity perturbations was eliminated by varying the 
. perturbation size until no changes in the coefficients occur for large 
. changes in perturbation size. Thus the half Sommerfeld condition, with 
both bloating-floating and fixed-fixed cases, can be examined. Principal 
stiffnesses are plotted vs. Sommerfeld number in Fig. 8 and cross-coupled 
stiffnesses in Fig. 9 where eccentricities are indicated at the top of 
se the graph. The half Sommerfeld floating-floating results agree well with 
the Reynolds fixed-fixed case. Figs. 10 and 11 give the damping co- 
efficients. The principal damping coefficients half Sommerfeld floating- 
floating curves again are very close to the Lund data points. For the 
cross-coupled damping, the Reynolds condition yields equal values of Buy 
and ae while the half Sommerfeld floating-floating curves are slightly 
above and below. Fig. 12 shows the stability threshold for the plain 


journal bearing. 


69 


tna Hee. eR 


ECCENTRICITY RATIO - © 


0.95 0.9 0.8 0.7 0.6 0.5 0.4 0.3 0.2 0.1 
100 


- Full Journal Bearing 
- L/D = 0.5 


2 Teg pT Ye ee OES 


rage: 


ue Half Sommerfeld Cavitation 


aoe §=Float-Float 


——e ao Fix-Fix 


e 
; Reynolds Cavitation 
oO Fix-Fix (Lund) 
i. 
< 10 


1.0 


DIMENSIONLESS STIFFNESS COEFFICIENT 


SOMMERFELD NUMBER 


i gives o% Princip:' Stiffnesses vs. Sommerfeld Number for Plain Journal Bearing 


70 


ECCENTRICITY RATIO - © 


0.8 0.7 0.6 O.5 0.4 


Full Journal Bearing 
L/D = 0.5 


Half Sommerfeld Cavitation 


—==> Float-Float 


—orm Fixed-Fixed 


adesynolds Cavitation 
O Fix-Fix (Lund) 


K_ (Negative) 
yx (Neg 


xy 
(Negative) 


e 
A 
isa) 
= 
oO 
_— 
1, 
tL, 
1] 
o 
Oo 
WY 
Ww 
wo 
Zz 
se 
bu 
al 
e 
n 
Y 
n 
A) 
4 
Zz 
oO 
eed 
n 
a 
a 
= 
— 
a 


SOMMERFELD NUMBER 


Figure 9 Cross-Coupled Stiffnesses vs. Sommerfeld Number for Plain 
Journal Bearing 
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Figure 10 Principal Damping vs. Sommerfeld Number for Plain Journal Bearing 
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Figure 11 Cross-Coupled Damping vs. Sommerfeld Number for Plain Journal Bearing 
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Figure 12 Stability Threshold vs. Sommerfeld Number for Plain Journal Bearing 
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In each case the half Sommerfeld floating-floating curves are very good 


agreement with the Reynolds fixed-fixed results. Only the cross-coupled 


a 


damping coefficients are somewhat different but generally by less than 

10%. These damping coefficients are of minimal importance in unbalance 
response of rotor-bearing systems and the effect on bearing stability, shown 

in Fig. 12, is also very small. The half Sommerfeld fixed-fixed values are 

not in good agreement with Lund's results. Thus, the half Sommerfeld floating~ 
floating condition is used to obtain the dynamic coefficients in this report. 
The half Sommerfeld condition avoids the iteration on the boundary loca- 


tion and saves a great deal of computer time. 


5.3 Reynolds Condition 

Reynolds condition with both fixed-fixed and floating~floating boundary 
locations were calculated with the series solution [Li, Allaire, Barrett, 
1977]. The only difference was the Reynolds boundary condition. Figs. 
13 through 17 gives the results for stiffnesses, damping, and stability. 
For comparison, the Lund perturbation results [Lund, 1969] are given 


as well, 


It is seen from the curves that the difference between fixed-fixed and 
floating-floating boundary locations is not nearly as large for Reynolds 
condition as it is for the half Sommerfeld condition. Some large differences 
occur in K_ at high Sommerfeld number where the value actually becomes 
negative for the floating-floating case. The cross-coupled damping B 

is also somewhat different at large S. The stability threshold actually 


has the largest difference. 
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Figure 13 Principal Stiffnesses vs. Sommerfeld Number for Plain Journal Bearing 
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Figure 14 Cross-Coupled Stiffnesses vs. Sommerfeld Number for Plain Journal 
Bearing 
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Figure 15 Principal Damping vs. Sommerfeld Number for Plain Journal Bearing 
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Figure 17 Stability Threshold For Plain Journal Bearing 
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5.4 Conclusions 


This section has presented numerous comparisons for both half Sommerfeld 

and Reynolds boundary conditions. The results indicate that the half 
Sommerfeld boundary condition differs from the Reynolds boundary condition 

by only a few percent in most cases (when the half Sommerfeld condition 

is properly applied to the dynamic coefficients). Thus the half Sommerfeld 
boundary condition has been used in this study. Although only the plain 
journal bearing has been examined in this section, it is believed that these 
results apply to pad bearings as well. On most pau. in a given bearing, the 
pad will have all positive pressures (due to a fully converging film thick-~ 

E% ness) or be all cavitated. Usually only one pad would even have a converging- 


diverging film thickness where Reynolds condition would be applied. It is thus 


se 


believed that the half Sommerfeld condition is adequate and much faster 


(in terms cf computer time) than the Reynolds condition. 
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